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A method is presented for computing approximate expressions for the stochastic 
force term £ 5 which appears in the Einstein-Langevin equation of stochastic gravity. 
Within this framework, £ a {, is a stochastic tensor field whose probability distribution 
mimics the probability distribution of the fluctuations of the quantum stress tensor 
operator; it is defined to be a random tensor field of zero mean whose correlation 
function is given by the expectation value of the symmetrized two point function of 
the stress energy fluctuation operator, called the noise kernel. Approximate expres- 
sions are obtained by means of a truncated Karhunen-Loeve transform defined on a 
random lattice of spacetime points. Due to the singular nature of the noise kernel, 
a coarse graining procedure is used to regulate divergences; as a result, the expres- 
sions obtained for £ a f, approximate values which might be seen by a probe measuring 
fluctuations in the stress energy using a sampling profile of finite width. Two real- 
izations of in Minkowski spacetime for the conformally invariant quantum scalar 
field in the Minkowski vacuum state are presented. 

PACS numbers: 04.62+v 



I. INTRODUCTION 

Stochastic gravity pQ is a theory first developed in the '90s with the intent to extend 
semiclassical gravity to include the effects of the backreaction of quantum fluctuations of 
matter fields on the geometry of spacetime. The mechanism by which this is accomplished is 
the introduction of a stochastic tensor field t; a b which, although classical in nature, neverthe- 
less replicates the probabilistic behavior of the quantum stress tensor fluctuation operator 
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tab = T a b ~ (Tab)- 

This fluctuation term is incorporated into the field equations for gravity via the Einstein- 
Langevin equation J2j 13] , 

G ab (x) = 8nG((f ab (x)) + U(x)) , (1.1) 

wherein the stochastic tensor field £ ao plays the role of the source driving the metric pertur- 
bations. Solutions to this equation will be a stochastic distribution of geometries determined 
by the properties of 

Up to its second moment, £ ab may be determined by the expectation value and two point 
correlation function of t a b, 

(U{x))s = (iab(x)) = (1.2a) 
(U(x)ZMx'))s = ({iab(x),i c 'd>(x')}) . (1.2b) 

If the fluctuations are assumed to follow a Gaussian probability distribution this is suf- 
ficient to uniquely determine ^ ao [2J. Thus, the critical object for which solutions must 
be found is the symmetrized two-point function of the stress tensor fluctuation operator, 
({iab(x), i c 'd'(x')}), which will hereafter be referred to as the noise kernel N a b c >d>(x, %')■ 

Since about 2000, the major focus of investigation in stochastic gravity has been on 
finding ways to compute usable expressions for the noise kernel in various spacetimes [IHH] . 
However, it is difficult to see how to extract physics directly from knowledge of the noise 
kernel. This paper is a first step in extending the theory, by deriving a method to compute 
explicit expressions for the stochastic tensor £ a & from the noise kernel itself. 

Unfortunately, due to the highly singular nature of the noise kernel, without some method 
of regularization an explicit functional form for £ ao does not exist; the stress tensor fluctu- 
ations will have infinite amplitude and be discontinuous at every point. This behavior is 
not expected to be physical, but is rather an artifact of the UV behavior of the field the- 
ory. To rectify this issue, we consider fluctuations not in the continuum limit but rather 
those which might be seen by some measurement apparatus of finite sampling width. This 
is approximated by coarse graining through the use of smearing functions, along the lines 
of [15]. 
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Explicit expressions for the coarse grained fluctuations are obtained by expanding the 
stochastic force (, a b via a Karhunen-Loeve transformation [IB], wherein £ ab is decomposed 
into a series of orthogonal basis functions with random, uncorrelated coefficients of known 
probability distribution, 

u*) = X>*2(*)- (1.3) 

i 

For the case of Gaussian fluctuations, each Z{ will be an independent, Gaussian random 
variable of zero mean, and the functions 4>ab( x ) are computed from the noise kernel. 



Typically, a series expansion along the lines of Eq. (1.3) will contain a countably infi- 
nite number of terms (one for each point in a countable, dense subset of the background 
spacetime), so in practice we find approximate expressions by truncating the series. The 
advantage of the Karhunen-Loeve transformation is that the series so generated minimizes 
the mean squared error introduced by the truncation procedure. As a general result, when 
£ a b is assumed to be Gaussian, each of the functions 0^ (x) may be approximated by a finite 
sum of functions which go like the noise kernel with one of the points fixed. 

As an unexpected advantage, this coarse graining procedure provides a natural means for 
using fluctuations to test the validity of semiclassical gravity. Although smearing functions 
introduce an artificial length scale into the theory, there is strong dependence on this length 
scale of the amplitude of the observed stress energy fluctuations. Correspondingly, for 
lengths scales in which the amplitude of these fluctuations exceed some threshold value with 
high probability one would expect the semiclassical approximation to no longer hold [TTHT9] . 

This paper is organized as follows. In Sec. [II], a brief review is given covering the relevant 
details of the Karhunen-Loeve transformation. In Sec. |III[ a computation of the smeared 



noise kernel in Minkowski spacetime is presented for the conformally invariant scalar field 



in the Minkowski vacuum state. In Sec. |IV[ we derive an expression for the coarse grained 
stochastic force in terms of the noise kernel, and present a modified Einstein-Langevin equa- 
tion which is potentially solvable. For the case of the Minkowski vacuum, two realizations 
of the stochastic force are computed numerically and displayed. Section [V] contains a sum- 
mary and discussion of the main results. Throughout this paper, units are used such that 
h = c = G = kB = l and the sign convention is that of Misner, Thorne, and Wheeler [20j . 
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II. KARHUNEN-LOEVE EXPANSION 

A random field R(x) with zero mean and correlation function C(x,y) may be expanded 
by means of a Karhunen-Loeve transformation [TB] as 



R(x) = Y,Z l r {i \x) , (2.1) 

i 

where the Z^s are centered, uncorrelated random variables and each r^(x) is a known 
(non-random) function of x. The functions r^'(x) are determined by the solutions to the 
eigenfunction equation 

J dyC{x,y)r®{y)=\ i r®(x). (2.2) 
The Zi's are determined, up to their second moment, by 

(Zt) a = (2.3a) 

(Z i Z j ) 8 = SijXi, (2.3b) 

where (...) s indicates that we are taking the stochastic expectation value (as opposed to the 
expectation value of a quantum operator). 



In general, the eigenfunction equation (2.2) is difficult to solve exactly and there will be 



an infinite number of eigenfunctions. However, it is possible to numerically approximate 



R(x) by discretizing to a lattice of points. Then the eigenfunction equation (2.2) becomes a 
much simpler eigenvector problem 

J2C(x u y J )r^(y J ) = X k r^(x l ), (2.4) 

3 

which may be solved using standard linear algebra techniques. Approximate eigenfunctions 
may be recovered via 

r (k \x) « X^^CfaxdrWfa) . (2.5) 

i 

When the correlation function C(x,y) is finite, continuous, and positive definite the series 



expansion in Eq. (2.1) converges in the mean-squared sense; that is, 



lim((R(x)-R k (x)) 2 ) s = , (2.6) 

k— >-oo 
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where 



R k {x) = Y J Z i r^\x) . 



(2.7) 



i=i 



Furthermore, the expansion procedure presented here has the benefit that it minimizes the 
mean-squared error introduced due to trunctation of the series, although the accuracy of 
such a discretization will depend on the density of points in the lattice. 

Additionally, if the random field R(x) is Gaussian in nature, then each coefficient Zi is 
an independent Gaussian random variable and the correlation function C (x, y) is sufficient 
to exactly specify all of the moments of each coefficient. In this case, the summation in 



Eq. (2.1) converges exactly (not just in the mean-squared sense). If the random field is 



non-Gaussian, then C(x,y) is in general not sufficient to uniquely specify each Z iy and the 
higher moments must be computed from the higher moments of R(x). 



Examining Eq. (2.5), it is apparent that each approximate eigenfunction is constructed 



from a finite sum of functions which look like the correlation function C(x, y) with one point 
fixed. We will call such functions 



= C(x,Xi) 



(2.8) 



shape functions. 

By reformulating the eigenvector equation as a matrix equation, 



C(x 1 ,x 1 ) ■■■ C(x 1 ,x N ) 



yCixN^x) ■■■ C(x N ,x N ) j 



A; 



(2.9) 



it is clear that each column of the correlation matrix C corresponds to one shape function 
^W. Let A be the unitary matrix of eigenvectors r % such that 



A T ■ C ■ A = diag(X) 



(2.10) 



Written thusly, it is clear by inspection of Eq. (2.5) that each approximate eigenfunction 
r( k \x) is given by 



(*). 



x) . 



(2.11) 
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Finally, we may write an approximate expression for the random field in terms of the eigen- 
values, eigenvectors, shape functions, and random coefficients as 

R(x)^J2 Z k X k 1 Ak^ {l) (x) . (2.12) 

i,k 

This provides a computationally inexpensive way to approximately obtain realizations for 
any Gaussian random field provided that the correlation function is known and satisfies the 
finiteness, continuity, and positive definiteness criteria. 



III. SMEARED NOISE KERNEL 



A general expression for the noise kernel for the conformally invariant scalar field was 
derived in [6] in terms of four derivatives acting on the square of the Wightman function, 

Nabdd' = Re {K abddl + g ab K cld < + g c >d>K ah + g a b9dd'K} (3.1) 

with 1 

9K a bc>d' = 4 (G-jb G ; d' a + G-(j a G-d'b) + G ]C >d> G- a b + GG- a bc'd' 

— 2 (G-b G ]C ' a d> + G- a G ]C >bd' + G-d> G. a b c ' + G- C i G ]a bd') 

+2 (G- a G-b R c 'd' + G.(j G t d' R a b) 

— (G- a b Rc'd' + G.tfd' R a b) G + -R c 'd' RabG 2 (3.2a) 



36.fC& _ 8 ^— G- p /b G- p a + G-b G- p r a p +G- a G- p ib p 

4 (g P G- a b p i — G- p i P G- a b ~ GG- a bpi P J 

—2 R' (2 G- a G : b — G G-ab) 

-2 (G. pl G. p ' - 2 G Rab - R' R ab G 2 (3.2b) 



p ,g 
r ;p g 



36K = 2G,p, q G p ' q + 4 (G- p > p ' G- q q + GG- 
-4 (G. p G. q ^' +G/G ;g v) 
+RG, pl G p ' + R'G )P G' P 

-2 (r G,/ + R' G/) G + ^RR'G 2 . (3.2c) 



1 Note that the superscript + on G + has been omitted for notational simplicity. 
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Primes on indices denote tensor indices at the point x' and unprimed ones denote indices at 
the point x. R ab and R c >d' are the Ricci tensor evaluated at the points x and x', respectively; 
R and R' are the scalar curvature evaluated at x and x' . 

Evaluated in Minkowski space, the noise kernel for this field in the Minkowski vacuum 
may be written as [21J, 



, T / /\ O a O h O c iG d i CTr a r} b )r c >Cr d ,) ^T}a(&Vd')b - VabVc'd' n , 

J^abc'd' [X, X) = -. — 2 —, — -. — -. . 3.3 

1 ' ; 48ttV 6 24tt% 5 1927t 4 (t 4 v ; 

Here (...) indicates symmetrization of the indices, 7] a b is the Minkowski metric, and rj ac ' = 
diag(— 1, 1, 1, 1) is the bivector of parallel transport in Minkowski space for Cartesian coor- 
dinates, a is the Synge world function, which in Minkowski space evaluates to 

a(x,x') = ^[-(t-ff + ix-x 1 ) 2 } , (3.4) 

and a a indicates the covariant derivative of a with respect to a. 



As can be seen by inspection of Eq. (3.3), the noise kernel diverges in the coincident limit 
and for null separations of points as l/(x — x') 8 . As a result of this divergence, no functional 



expression for £ a ;, exists; each eigenvalue in Eq. (2.2) will be infinite. Similar to the case of 
forces in Brownian motion, £, a b{x) will be infinite and discontinuous at every point. 

In order to rectify this problem, it is necessary to find a method to regulate the noise 
kernel. The method used here is coarse graining through the use of smeared operators, 
following the procedure outlined in [T3]. For our purposes, it is sufficient to smear only 
along the time coordinate. 2 

Using a Lorentzian smearing function 

wm = wh^) • < 3 - 5 > 

we compute the smeared Wightman function in Cartesian coordinates as 
G(x, x') = J dt" df G(x, f, t, t')W(t" - t)W(f - f) 



At 2 - 4« 2 - r 



2 



47r 2[ r 4 _ 2r 2 (At 2 - 4a 2 ) + (At 2 + 4a 2 ) 



(3.6) 



2 Spatial-only smearing is also sufficient to regulated the divergences, as is shown in [3]. 




FIG. 1: -/Voooo component of the unsmeared noise kernel in the x-t plane (y = z = 0). 

where At = (t — t') and r = y(x — x') 2 . The resulting expression for the iVoooo component 
of the smeared noise kernel is 3 



N 0000 (x,x') 



x 



127r 2 [r 4 - 2r 2 (At 2 - 4a 2 ) + (At 2 + 4a 2 ) 2 ] 6 
"(r 2 - At 2 ) 6 (3r 4 + 10r 2 At 2 + 3At 4 ) 

-32a 2 At 2 (r 2 - At 2 ) 4 (33r 4 + 42r 2 At 2 + 5At 4 ) 

-64a 4 (r 2 - At 2 ) 2 (9r 8 + 132r 6 At 2 - 526r 4 At 4 - 212r 2 At 6 + 21 At 8 ) 

-512a 6 (6r 10 - 95r 8 At 2 - 552r 6 At 4 + 390r 4 At 6 - 246r 2 At 8 - 15At 10 ) 

+512a 8 (45r 8 + 1300r 6 At 2 + 1326r 4 At 4 + 436r 2 At 6 + 221At 8 ) 

+8192a 10 (36r 6 + 309r 4 At 2 + 110r 2 At 4 + 57At 6 ) 

+49152a 12 (25r 4 + 78r 2 At 2 + 21At 4 ) 

+131072a 14 (I8r 2 + 13At 2 ) 

+ 1769472a 16 ] . (3.7) 



In the a = limit this reduces to 



N, 



0000 



3r 4 + 10r 2At 2 + 3At 4 
127r 2 [r 2 - At 2 ] 6 



{3.1 



which agrees exactly with Eq. (3.3) 



Figures [T] and [2] show the comparison between the N 0000 component of the unsmeared 
and the smeared noise kernel; in the smeared case all divergences have vanished and the 



3 For brevity, only the iVoooo is displayed. The expressions for the other components are similar, but 
depending on the component may have explicit dependence on separations in the x, y, and z coordinates. 
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0.5 1.0 



FIG. 2: iVoooo component of the smeared noise kernel in the x-t plane (y = z = 0) for a = 0.25. 



function is continuous and finite everywhere. 



IV. GAUSSIAN STRESS ENERGY FLUCTUATIONS 



The goal is to find a Karhunen-Loeve expansion for the stochastic tensor field £, a bi x ) 



of the form given in Eq. (1.3) following the scheme outlined in Sec. II In this case, the 



eigenfunctions 0„ 6 s are given by solutions of the eigenfunction equation 



AV-sii'liv,,"!!,/) 



ab \ X ) > 



(4.1) 



where the only change from Eq. (2.2) is the addition of the sum over the primed indices. As 



before, if the fluctuations are assumed to have a Gaussian probability distribution, then the 
noise kernel is sufficient to uniquely specify £ a 6. 

Discretizing to a lattice of points gives the eigenvector equation 

N 

£ N **^ ^A x j) = A* <&\*i) . ( 4 -2) 

3=0 

which can be solved using standard linear algebra techniques. In order to prevent any 
artificial periodic structure from appearing as a result of this discretization procedure, the 
points in the lattice are chosen at random from a closed volume of spacetime. As the error 
in the approximation scheme is related to the density of points, it is therefore possible that 
regions of higher error will occur due to a locally underdense sampling of points; however, 
this issue can be controlled by the choice of a sufficiently large number of points. 
As in Sec. [IT], shape functions are constructed in terms of the noise kernel, giving 



(0 AT c'd' 
C c'd' IX ab 



X, Xj) . 



(4.3) 



10 



However, in contrast with Eq. (2.8), here it is necessary to introduce constant tensors to 



pick out individual components of the noise kernel. For simplicity, choose each such that 
a single component is equal to one and all other components are zero. For instance, choose 
c^J = 5^5%. Thus, there are 16 functions ^ a b(x) for each point X{ in the lattice; however, due 
to the symmetries present in the noise kernel only nine of these are linearly independent. 
Reformulating: clS du matrix equation gives 



Nj d ' 



N ab cd (x u x N ) 



\ 



1 4>%{*i) ^ 



\ N ab C ' d '( X NiXl) ■■■ N ab c ' d '(x Nl X N ) J \^1,(X N ) J 

1 ^ 



(4.4) 



where each N, 



c'd' 



ab 



X 



,Xj) represents a 16 x 16 block of the matrix (and each 



16 element block of the vector). As before, each column of the above matrix for the noise 
kernel represents one of the shape functions (x) evaluated on the points of the lattice. 
The final expression for ^ a b{x) becomes 4 



t,ab{x) 



Y,Z k K l A km+l) ^{x), 

j,k,l 



(4.5) 



where once again we define A to be the matrix of eigenvectors of N abc id' ■ 

We are now in a position to write down the smeared Einstein-Langevin equation, 



Gah(x) 



8ttG 



(f ab (x)) +^Z fe A^U fci(16i+0 *S' ) (x) 



j,k,l 



(4.6) 



Solutions to this equation will involve a stochastic distribution of geometries; however, 
within the linearized gravity approximation, for small fluctuations (large a) it may be suf- 
ficient to compute the metric perturbation generated by each shape function ^^'^ individ- 



ually. Schematically, one may solve Eq. (4.6) for the metric perturbations induced by the 



4 The purpose of the 16j + I index of eigenvector matrix A is to distinguish between the 16 possible shape 
functions associated with each point of the lattice. 
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FIG. 3: One realization of £ a b along the x-t plane (y = z = 0) with a = 0.25 for a random lattice 
of 50 points chosen within the hypercube bounded by [—0.5,0.5]. 

stochastic force using [18] 



where G ah c ' d ' (x , x') is the semiclassical retarded Green function for the operator acting on 



metric pertubations. For maximally symmetric spacetimes such as Minkowski space, one 
would need to compute backreactions for only nine unique functions, due to the translation 
symmetry of the shape functions. 

Even without solving the backreaction problem, it is possible to discern some of the 
behavior of the fluctuations by examining realizations of the stochastic force £ a fc. The most 
immediately apparent of these is the strong dependence of both the amplitude and the 
frequency of the fluctuations on the smearing parameter a. As shown in Figs. [3] and [4| as 
a decreases the amplitude and frequency of the fluctuations increases, until in the a = 
limit both the amplitude and frequency of the fluctuations are infinite. This is consistent 
with forces in the case of Brownian motion, where in the short distance limit the forces are 
infinitely strong and discontinuous everywhere. 

As a result of this dependence of the amplitude and frequency on the smearing parameter 
a, the smeared noise kernel provides a natural length scale for use in evaluating the validity 
of the semiclassical approximation; since the semiclassical approximate is suspected to not 
be valid when fluctuations of the stress energy tensor are large, then it is unlikely to be 




j,k,i 



(4.7) 



h a b. Thus, the set of shape functions vErC^O may be used to generate a basis for the induced 
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FIG. 4: One realization of £ a ft along the x-t plane (y = z = 0) with a = 0.1 for a random lattice of 
50 points chosen within the hypercube bounded by [—0.5,0.5]. 

valid on length scales smaller than some «o chosen such that when a > «o the fluctuations 
remain below a given threshold value with high probability. 

A second feature evident from Figs. [3] and [4] is that fluctuations tend to propagate along 
null vectors, as would be expected from the fact that the noise kernel is large for points 
near null separation; however, the length scale of this propagation is dependent on a, with 
small a leading to shorter propagation lengths due to the higher frequency of interfering 
fluctuations. 

V. DISCUSSION 

In this paper, we have derived a method for generating approximate, coarse grained 
expressions for the stochastic force term ^ which appears in the Einstein-Langevin equation 
of stochastic gravity. This method makes use of a Karhunen-Loeve transform to express 
£ a b in terms of decorrelation random coefficients and known non-random functions of the 
spacetime points. In order to render this method amenable to numerical computation, the 
Karhunen-Loeve transform was discretized to a random lattice of points. 

Coarse graining proved necessary due to the highly singular nature of the noise kernel 
and the resulting pathological UV behavior of the fluctuations. While this behavior is 
not expected to be physical, it does demand that some method be chosen to regulate the 
noise kernel. The method chosen here was smearing along the time direction by use of a 
Lorentzian smearing function with an associated length parameter a; this approximates the 
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stress energy fluctuations which might be seen by a probe with a sampling profile which 
matches that of the smearing function. 

The expressions for the stochastic force generated by this procedure are composed of 
finite summations of shape functions which look like the coarse grained noise kernel with 
one point fixed multiplied by random coefficients. When the fluctuations are assumed to be 
Gaussian in nature, each of these coefficients is also Gaussian. 

As a result of the dependence of the amplitude of the stochastic force upon the smearing 
parameter a, this method provides a natural means by which the validity of semiclassical 
gravity may be tested. Since semiclassical gravity is expected to fail when fluctuations are 
large, then it can be assume to be valid only for length scales greater than some ao chosen 
such that when a > ao the fluctuations remain below a given threshold value with high 
probability. 

Several realizations of £ a & have been computed in Minkowski spacetime for the conformally 
invariant scalar field in the Minkowski vacuum state, and two such realizations have been 
displayed. By investigating these realizations it is possible to discern some of the behavior 
of the stress energy fluctuations. For the case of the conformally invariant scalar field, one 
sees that fluctuations tend to propogate along null vectors, as one would expected from the 
fact that correlations in the stress energy for that field are large for nearly-null separations 
of points. However, this propogation length appears to be tied to the smearing parameter 
a due to the higher frequency of interfering fluctuations. 

Finally, the method presented here provides a pathway to finding explicit solutions of the 
coarse-grained Einstein-Langevin equation. Such solutions would be a stochastic distribution 
of geometries whose properties are determined by the shape functions and the random 
coefficients Z{. The author hopes to revisit this problem and compute such solutions in 
future work. 
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